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Abstract. We prove the existence of the hmiting spectral distribution (LSD) of 
I symmetric triangular patterned matrices and also establish the joint convergence 

of sequences of such matrices. For the particular case of the symmetric triangular 
Wigner matrix, we derive expression for the moments of the LSD using properties 
^ of Catalan words. The problem of deriving explicit formulae for the moments of 

the LSD docs not seem to be easy to solve for other patterned matrices. The LSD 
of the non-symmetric triangular Wigner matrix also does not seem to be easy to 
establish. 



^ 1. Introduction 

^ Two of the most important problems in operator theory are the invariant subspace 

I ^1 problem and the hyper invariant subspace problem. Let "H be a separable Hilbert 

^ space (infinite dimensional) and let i3('H) be the algebra of bounded operators on Ti. 

^ An invariant subspace of A e B{H) is a subspace C Ti such that A('Ho) C Tio 

CN and a hyper invariant subspace of A is a subspace Hq C H that is invariant for every 

^ operator B G i3(H) that commutes with A. The invariant subspace problem asks 

T-H whether every operator in Bl^-L) has a closed non-trivial invariant subspace. The 

O hyper invariant subspace conjecture states that every operator in i3('H) that is not 

^ a scalar multiple of the identity operator has a closed, non-trivial hyper invariant 

y—i subspace. For many years there were attempts to prove this result. Recently there 

j>! have been attempts to disprove this conjecture. 

A natural candidate which stood out for counter examples was the DT-operators 
^ defined by Dykema and Haagerup [11]. This is related to the following asymmetric 

^ matrix defined and studied in their work. It is a triangular version of the celebrated 

Wigner matrix (which is a symmetric matrix and, in its simplest form, has i.i.d. 
entries). 



Date: 5th October, 2011. First revision January 26, 2012. Second revision April 06, 2012. 

2000 Mathematics Subject Classification. Primary 60B20; Secondary 46L54. 

Key words and phrases. Triangular matrices, Wigner, Hankel, Toeplitz and Symmetric Circulant 
matrices, limiting spectral distribution, asymptotically free, DT-operators, Catalan words, symmet- 
ric words, semicircular law, joint convergence of matrices. 

R. Basu's research is supported by Loeve Fellowship, Department of Statistics, University of 
California, Berkeley. 

A. Bose's research supported by J.C.Bose National Fellowship, Dept. of Science and Technology, 
Govt, of India. 

1 



2 



R. BASU, A.BOSE, S. GANGULY, AND R. S. HAZRA 







^1,2 


tl,3 • 




tl,n 







^2,2 


^2,3 • 


■ i2,n-l 


^2,ra 


T = 








^3,3 • 


■ i^3,n-l 


^3,71 










. 








where {tij)i<i<j<n are i.i.d. complex Gaussian random variables having mean and 
variance 1/n. 

Let /X be a compactly supported measure on the complex plane C and let Dn be 
a diagonal matrix with /^-distributed independent random variables which are also 
independent of T„. Let r„ denote the functional ^ Tr. Dykema and Haagerup [H] 
showed that Z„ = cn~^/^T„ + Dn converges in *-moments, that is, 



[1.2) 



exists for every A; G N and for all e{l) , e{2) , . . . , e{k) G {*,!}• Further, the lim- 



its DT{n,c) in (L2) are operators and may be viewed as elements of a W* non- 
commutative probability space (A^,r) where is a von Neumann algebra and r 
is a normal state. They also showed that these DT-operators are decomposable (an 
operator K is decomposable if for every cover C = f/ U F of the complex plane of 
open subsets U and V there are i^-invariant closed subspaces H' and H" such that 
a{K\y^i) C U and o"(i^'|-^") C V and Ti = Ti' + Ti"). As a consequence, it was shown 
in Dykema and Haagerup [HI |12j that DT-operators whose spectra contains more 
than one point have non-trivial closed hyper invariant subspaces. Thus the study of 
the spectrum of the operator DT is of immense importance. 

Some properties of the *-limit of n~^^'^Tn were derived by Dykema and Haagerup 
[TT| m\ . Sniady [TH] and Shlyakhtenko \L7\. In particular, considerable analytical 
techniques as well as combinatorics of operator valued free probability theory were 
used to show that 



1 

E^Tr ( ^ 
n \ n 
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Moreover, the sequence {pr^rxyy} are moments of a probability measure v supported 
on [0, e] and given by 



dvi^x) = ip{x)dx where ip : (0, e) — t- 



:l4) 



is the unique solution of 



f sinv . 
ip expffcotf 



TT 



sint>exp(— f cot v). 



:l5) 



The plot of the above function is given in figure [TJ 
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Figure 1. Plot of function V' 



Now consider the triangular symmetric Wigner matrix: 



n 



Xl2 



Xl2 
X22 



Xl3 
X23 



■^1(71— 1) -^In 
X2{n-1) 



a;i(n-l) X2{n-1) 
Xln 













(1.6) 



where {xij} are independent and identically distributed random variables. 

Note that if we remove the symmetry from the above matrix and consider the 
matrix with i.i.d. entries above the anti-diagonal, then it can easily be linked to T„. 
If we multiply this matrix by a matrix P whose entries on the anti-diagonal are 1 
and the rest are zero, then we get a matrix of the form T„. It can also be seen that 
{W^y = TnT* where T„ is an upper triangular matrix (1.1) but has dependent rows 



and columns now. 

Our goal is to study the spectral properties of n'^^'^W^. More generally, we consider 
symmetric triangular patterned random matrices and whose entries are not necessarily 
Gaussian. 

We use the classical moment method to show that the limiting spectral distribution 
(LSD) of triangular Wigner, Toeplitz, Hankel and Symmetric Circulant matrices exist 
when the entries have mean zero and variance one and, are either (i) independent and 
uniformly bounded or (ii) i.i.d. 
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The joint convergence of p sequences of random matrices is defined by the con- 
vergence of ^ETr(g) for all monomials q of the p matrices of order n. The LSD 
result for one patterned matrix at a time can be extended to the joint convergence of 
any combination of these patterned matrices when we assume that the sequences are 
independent and the entries are independent with mean zero and variance one and 
have uniformly bounded moments of all orders. All the limit results are universal (do 
not depend on the underlying distribution of the entries). 

The identification of the LSD of general patterned matrices is not a trivial problem. 
For certain (full) circulant type matrices this can be done (see for example Bose and 
Mitra [4J, Kologlu et al. [I4j, Massey et al. [15]). One may recall that for the full 
Toeplitz and Hankel matrices, no formulae of any sort are known for the LSDs or 
their moments (see Bose and Sen [5], Bryc et al. [9], Hammond and Miller [13]). 
By exploiting the Catalan recursions inbuilt in the Wigner matrix, we show that 
the (2/c)th moment of the LSD for is given by (1.3) and the odd moments are 
zero. However, moment formulae for other matrices remain elusive and appear rather 
difficult to obtain. 

Identification of any structure in the joint limit moments appears much harder, 
even for the triangular Wigner matrix. Unlike the full Wigner, where independent 
sequences of Wigner matrices are asymptotically free, the triangular Wigner matrices 
are not so. It would be interesting to unearth the structure of the joint limit for the 
triangular Wigner matrices. Likewise, the full Reverse Circulant joint limit satisfies 
half independence with symmetric Rayleigh marginals and the full Symmetric Circu- 
lant joint limit is totally independent with Gaussian marginals (see Bose et al. [S]). 
However, the joint moment structure for their triangular versions seems rather hard 
to obtain. We do not even know the LSD or the moments of the marginals. 

In Section[2]we describe our setup, notation and state and prove the main existence 
results for the LSD. In Section |3] we concentrate on the triangular Wigner matrix and 
derive the moments of the LSD. In Section |4] we outline how our results may be 
extended to joint convergence and also discuss many interesting questions which arise 
in the study of random triangular matrices. 



2. Patterned triangular matrices and method of moments 

For recent developments on random patterned matrices we refer the readers to Bose 
et al. [7]. We adopt the method developed in Bryc et al. |9j and Bose and Sen 0. 
Here is a quick description of the concepts and notation we borrow from these works. 

A sequence or bi-sequence of variables {xi^i > 0} or {xif, i,j > 1} will be called 
an input sequence. Let Z be the set of all integers and let Z+ denote the set of all 
non-negative integers. Let 

L„ : {1, 2, . . . ny ^Z'^, n>l, d = l,2 (2.1) 

be a sequence of functions such that L„+i(z, j) = Ln{i,j) whenever I < i, j < n. We 
shall write Ln = L and call it the link function and by abuse of notation we write 
as the common domain of {Ln}- Patterned matrices are matrices of the form 



Xn = {{xL(ij))). 



(2.2) 
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Throughout this article we assume that the triangular version is derived from a pat- 
terned random matrix X^. Let -'^"-be the upper triangular version of the matrix 

where the (i, j)-th entry of X" is given by xm^j) if (i+j) < n + 1 and otherwise. We 
shall often drop the subscript n, and simply denote the matrix by X". The triangular 
Wigner matrix has already been defined. Here are examples of the triangular versions 
of some more common patterned matrices: 



1. Triangulcir Hankel matrix. 



X2 
X3 



X3 X4 
X4 



Xn -^n+l 
Xn+1 



Xn Xn+l 
Xn+1 


















2. Triangular Toeplitz matrtix. 



'J-'U 



Xo 


Xl 


X2 . 


Xn—2 Xn—1 


Xl 


Xo 


Xl . 


■ Xn-3 


X2 


Xl 


Xo . 


. 



x„_i ... 
3. Triangulcir Symmetric Circulant Matrix. 



QU 



Xo Xl X2 

Xl Xq Xl 

X2 Xl Xo 





X2 Xl 

xs 







(2.3) 



(2.4) 



(2.5) 



It is to be noted that a triangular Reverse Circulant matrix is the same as a triangular 
Hankel matrix. 

To prove the existence of LSD via moment method, it suffices to verify the following 
three conditions. 



(1) For every k > 1, 

(2) For every k > 1, 



— >■ as n — >■ 00 (Condition Ml). 



E 



iT.((-l=X:rt-E(lT.-((-lx„«)- 

(3) Ylk^2k^^'' — °° (Carleman's Condition). 
The two main conditions we need are: 



O [n ^) as n ^ 00 (Condition M4). 



Assumption A: {xi} or {xij} are independent with mean and variance 1 and either 
(i) have uniformly bounded moments of all orders or (ii) are identically distributed. 
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The restriction of the so called Property B of Bose and Sen [5] to the triangular 
matrices is given below. 

Assumption B: The link function L satisfies Property B', i.e., 

A = sup sup sup #{/ :l</<n, A; + /<n + l, L{k, I) = t} < oo. 

n t l<k<n 

Since the moment method entails computing all moments, we shall also use the fol- 
lowing assumption. 

Assumption C: {xi} or {xij} are independent with mean 0, variance 1, and are 
uniformly bounded. 

In the course of verifying conditions (1), (2) and (3), we shall need the notion of 
circuits, words, etcetera from Bose and Sen [5]. 

Any function vr : {0,1,2, h} — )• {1,2, ■■■ ,n} with 7r(0) = 7r(/i) is called a 
circuit of length h. The dependence of a circuit on h and n will be suppressed. 

Two circuits tti and 712 are equivalent if and only if their L values respectively 
match at the same locations. That is, 

{L{n,{i - 1), 7ri(z)) = L{n,{j - ^ L{7iS - 1), ix^it)) = Hrt^U - 1), 7r(j))}. 

Any equivalence class can be indexed by a partition of {1,2, ■ ■ ■ , h}. Each block of 
a given partition identifies the positions where the L-matches take place. We can 
label these partitions by words of letters where the first occurrence of each letter is 
in alphabetical order. For example ii h = 5 then the partition {{1, 3, 5}, {2, 4}} is 
represented by the word ababa. This identifies the circuits vr such that L(7r(0), 7r(l)) = 
L(7r(2),7r(3)) = L(7r(4), 7r(5)) and L(7r(l), 7r(2)) = L(7r(3), 7r(4)). Let w\i] denote the 
i-th entry of w. A word is said to be pair-matched if each letter appears exactly twice. 
Denote by the length of the word w. 
For a pair-matched word w of length 2k, let 

U*Aw) = {iT:w[i]=w[j]^L{7r{z-l),7r{t)) = L{7r{j-l),7r{j))} 
Ux{w) = {vr : = ^ L(7r(z - l),7r(i)) = L(7r(j - l),7r(j))}. 

Note that Ilx{w) is nothing but the equivalence class of circuits corresponding to the 
word w. We shall write now 

1 ^ 

rPfc72 5Z Yl^L{iTii-l)Mi))'^Mi-'^}+-^ii)<n+l} 
TT-.n circuit i=l 

1 ^ 

5Z W^L{-n{i-l),-n{t))'^{7,{i-l)+iT{t)<n+l)- 
w 7ren{«)) i=l 

Let us denote 

k 

2;L(7r(i-l),7r(i))l{7r(j-l)+7r(i)<n+l}- 

i=l 



-Tr((^X 
n \ n 



u\k\ 
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For a fixed fc, let us define the class 11^ as follows. 

= {tt : TT is circuit of length k, 7r(z — 1) + 7r(z) < n + 1,1 < i < k}. 
Let us denote 

U,,xiw) = Uxiw)f]U-,, and Ul^w) = n^(^) fjn^,. 
For every pair-matched word w of length 2k, define, if the limit exists, 

Pu,x{w) = hm — — 

where for any set A, jj^A denotes the number of elements of A. 

When there is no chance for confusion, we simply write Pu{w) for Pu,x{w) and 
Hiiw) and ^\{w) for Ili xiw) and n^j^(w) respectively. 

By using the arguments of Bose and Sen [5], we immediately get the following 
theorem. 

Theorem 2.1. Let {X^} he a sequence of patterned triangular matrices satisfying 
Assumptions A and B. Suppose for every k > 1 and for every pair-matched word w 
of length 2k, Pu{w) exists. Then the LSD of ^ exists a.s. The LSD is universal, 
symmetric about and is determined by the even moments 



^2k= J2 k>l. (2.6) 

w pair-matched, 
\w\=k 

Proof. We provide a brief sketch of the proof. By using a truncation argument (see for 
example Bose and Sen 0), we may work under the stronger Assumption C. The first 
lemma in particular shows that the odd moments are zero and the LSD is symmetric. 

Lemma 2.1. Suppose Assumptions B and C hold. Let w be a matched word of length 
k which is not pair-matched. Then 

hm = 0. 

Proof. First note that if w is non- matched then E[X7r*] = since the entries have 
mean zero. Now consider a w which has more than two matches. Let = 
riiLi ^L(7r(i-i),7r(i))- Clearly then, < \Xt,\ and hence E|X^*| < E|X^|. It is 

proved in Bose and Sen [5] that 

hm J2 El^^l = 0. 

The lemma follows immediately from this since n^(w) C ir{w). □ 

The next lemma helps to verify Condition M4. The proof is a direct generalization 
of Lemma 2 of Bose and Sen [S]. We skip the details. 
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Lemma 2.2. Let {X^} be a sequence of random patterned triangular matrices satis- 
fying Assumptions B and C. Then 

1 4 



E 



-Tr(^X:)'=-E(-Tr(^X„")'=) 



O . 



Since there does not exist any pair-matched word of length for k odd, it follows 
from Lemma 2.1 that 

1 VU 

(3k = lim E(-Tr(-^)'=) = 0, if A; is odd. 

n-S>oo n y/n 

It also follows that 

(32k = lim E{-Ti{^y') = V p^{w). 
n-s>oo n Wn ^ — ^ 

w pair-matched, 
\w\=k 

Hence Condition Ml also holds. From Lemma 2.2 it follows that Condition M4 is 
satisfied. It remains to show that the sequence {Pk} satisfies Carleman's condition. 
It follows easily from Assumption B that for every pair-matched word of length 2k, 

Hence Pu{uj) < A*^ and so l32k < w^^'' ^^'^ Carleman's condition is thus easily 
verified. □ 



We can use Theorem |2.1| to prove the existence of LSD for the triangular matrices 
listed earlier. 

Theorem 2.2. Let {X^} be any of the following triangular matrices with input se- 
quence satisfying Assumption A: triangular Wigner, triangular Hankel, triangular 
Toeplitz and triangular Symmetric Circulant. Then the LSD for ^ exists almost 



surely. The LSDs are universal and are symmetric with even moments given by (2.6) 



Proof. As before, without loss, we shall work under Assumption C. Notice that the 
link functions of all these matrices satisfy Assumption B. Thus Carleman's condition 
and Condition M4 follow immediately if we can establish Condition Ml. The combi- 
natorics to show the existence of Pu{w) is done case- by-case. However, all the proofs 
are similar and use the volume method. We provide a sketch only for the triangular 
Wigner and the triangular Hankel matrices. In the full matrix version of the matrices 
considered here, p{w) is evaluated as an integral of an indicator function (see for ex- 
ample Bose et al. [7j). Here also, Pu{w) turns out to be an integral, but of a different 
indicator function corresponding to H^(u;) instead of Il*{w). 

Triangular Wigner Matrix: Let Lw denote the Wigner link function. Following 
Bose and Sen [5j, a word is said to be Catalan if it is pair- matched and deleting all 
double letters leads to the empty word. For example abba is Catalan and abab is not. 
Now it is shown in Bose and Sen [5] that if w is not Catalan then, 

p{w) = hm — — = 0. 
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As Il*{w) I) Ill{w), it follows that p„(w) = if ty is non-Catalan. Thus we now focus 
on Catalan words of length 2k. We shall call any 7r(i) (or by abuse of notation any 
i) a vertex. Any 7r(i) is a generating vertex if a letter appears for the first time at i 
(when read left to right). We also define 7r(0) (or 0) to be a generating vertex. 

Lemma 2.3. Let w be a Catalan word of length 2k. Let S denote the set of all 

generating vertices of w. Then for all j ^ S, there exists a unique i & S such that 
i < j and 7r(j) — 7r(i) for all tt e Il*{w). 

Proof. Let tt G Il*{w). Let j be the minimum index of a non-generating vertex of w. 
Clearly then, w[j — 1] = w[j] and hence 7r(j — 2) = vr(j). Since j > j — 2 E S , the 
result is true in this case. Now let j be any non-generating vertex. Let us assume that 
for every non-generating vertex with index less than j, the result holds. Let i < j 
be the index of first occurrence of w[j]. Let Wi be the subword formed by letters 
between w[i] and w[j]. Since w is Catalan, wi is also Catalan, and it can be easily 
shown that, 7i{i) = 7r(j — 1) and hence 7i{j) = 7r{i — 1). If (i — 1) G S, then we are 
already done. If (i — 1) ^ S, then also by induction hypothesis the result holds. 

If the i corresponding to a fixed j is not unique then we have a non-trivial re- 
lation between two generating vertices which implies #11* (w) = 0{n^), and hence 
contradicting the fact that lim„_^oo ^^i+r^ — 1- This proves the uniqueness. □ 

Definition 2.1. For any j (not necessarily in S), let us denote by 0(j) the unique 
vertex such that 



Note that if j G S. then j = 4>{j). Next we note that among the 2k equations, 
n{i — 1) + 7r(i) < n+1, 1 < i < 2k, each equation is repeated twice, as w[i] = w[j] =^ 
7i{i — 1) + 7r(i) = 7r(j — 1) + 7r(j). So we can write 



U.l{w) = {tt : w\i] = w[j] Lw{n{i - 1), 7r{i)) = Lw{TT{j - 1), 7r(j)), 
7r{i - 1) + 7r(i) <n + l\/ieS - {0}} 
= {tt : it{j) = 7r(0(j)) Vj ^ S, 7r(0(i - 1)) + 7r(0(z)) < n + 1 Vz G 5 - {0}}. 



#n^(w) = #{(-"0, ■■■,V2k) -VieUn VO < ? < 2k, Vi = Vi ^ S", 

^^<^(i-i) + v^(^i^ <l + l/n Wi e S - {0}, vq = V2k} 
= -.VieUn Vi G 5*, ^^(j-i) ^^(i) < 1 + 1/n Vi G 5* - {0}}. 

Prom the above equation it follows that is nothing but the Riemann sum for 

the function 



4>{j) G S, 4>{j) < j and 7r(j) = 7r(0(j)) for aU tt G U*{w). 



Now we use the standard volume method arguments. Let us define 




(2.7) 



Then, 



Iw{vs) = -^(^^</.(i-i) + 'v^pit) <l,ie S - {0}) 



10 R. BASU, A.BOSE, S. GANGULY, AND R. S. HAZRA 

over [0, 1]'^"'"^. Since the function is clearly Riemann integrable, the Riemann sum 
converges to the integral 

lim ^#n*(^i;) = / Iwivs)dvs. (2.8) 
It follows that Pu{w) = lim„_s.oo ^TTF#n*(?j7) exists. 

Triangular Hankel Matrix: The Hankel link function is L{i,j) = i + j. Here 

n^(t/7) = {n : w[i] = w[j] =^ 7i{i — 1) + 7r(i) = 7i{j — 1) + 7r(j), 7r(z — 1) + 7r(z) < n + 1}. 

Let S denote the set of all generating vertices of w. For every i G S* — {0}, let jj 
denote the index such that w[i] = w[ji]. Let us define Vi, Un,vs as in (2.7). Then 

#n*(w) = #{(t;o, ■■■,V2k) ■VieUn,y 0<i< 2k,V(^i_i) + Vi = + Vjyi G S, 

+ Vi < 1 + l/n,Wi e S - {0},vo = V2k}- 

It can easily be seen from the above equations (other than vq = V2k) that each of 
the {vi : i ^ S} can be written uniquely as an integer linear combination Li{vs)- 
Moreover, Li{vs) only contains generating vertices of index less than i with non-zero 
coefficients. For all i E S, let us define Li{vs) = fj. 
Clearly, 

#n*(w) = if{ivo,...,V2k)■■VieUnyO<i<2k,vo = V2k,v^ = Li{vs)y^^S, 

Vi-1 + Vi <l + l/n^i e S - {0}} . 

Integer linear combinations of elements of f/„ are again in f/„ iff they are between 
and 1. Hence, 

#Ul{w) = #{vs:v,e Un G 5, t;o = L2k{vs)^ < L\{vs) < 1 W ^ 5, 

Li^i{vs) + L,{vs) < 1 + 1/n G S - {0}} . ^^'^^ 



From (2.9) it follows that ^^i}^'' is nothing but the Riemann sum for the function 

Ih{vs) = 1(0 < L\{vs) <l,iiS,vo = L'^k{vs),Li^i{vs) + U{vs) < 1 G S - {0}) 

over [0, 1]*^+^. Since the function is clearly Riemann integrable, the Riemann sum 
converges to the integral 



lim -i-n:#nKw) = / lH{vs)dvs. 



Hence p„(w) exists for every pair-matched word w. 

Incidentally, following Bose and Sen [5] , a pair-matched word is said to be symmet- 
ric if every letter appears once in an even and once in an odd position. For a Hankel 
matrix, p{w) = if w is not a symmetric word. Hence it also follows that Pu{w) = 
for every such word w for the triangular Hankel matrix. □ 

The LSD of the Reverse Circulant, the Symmetric Circulant, the Toeplitz and the 
Hankel matrices are unbounded. The same is true for their triangular versions. 

Theorem 2.3. The LSD of triangular Hankel, Toeplitz and Symmetric Circulant 
matrices have unbounded support. 
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Proof. We shall prove the theorem only for the triangular Hankel matrices. The proof 
for other matrices is similar and is omitted. Recall that the 2/c-th moment of the LSD 
for the full and triangular versions are given by 



w pair-ri 
\w\ 

We claim that for every w, 



w) 



w pair-matched, w pair-matched, 

|u;|=fc |w)|=fc 



Pu{w) > ^p{w). 

To see this, first note that using the already introduced notation, 



Puiw)= / lHivs)dvs. (2.10) 

^[0,lf+i 

Also, clearly, by a similar logic for the full Hankel matrix, 

p{w) = [ J(0 < L\{vs) <l,t^S,Vo = L'2,{vs))dvs. 

Making a change of variable yi = Vi/2 for all i & S, since L- are linear functions, we 
see that 



2^ 



'[0,l/2]'=+i 

Now it is trivial to note that 



/ /(O < L\{ys) < 1/2,2 ^S,yo = L'^k{ys))dvs 

J[0,l/2]'=+i 

K{vs)dvs say. (2.11) 



K{vs) < Ih{vs) for all vs- 



The claim now follows from (2.10) and (2.11). As a consequence, 

P2k > 

and the result now follows from noting that the LSD of Hankel matrix has unbounded 
support. □ 



3. The LSD of Triangular Wigner Matrices 

As mentioned earlier, in general the LSD is hard to identify and we have the answer 
only for the triangular Wigner matrix. For the full Wigner matrix, p[w) = 1 for each 
Catalan word, and as a consequence the LSD is the semicircle law. In the following 
table we list the values of Pu{w) for Catalan words of small lengths. It may be ob- 
served that the contributions are equal within certain isomorphic classes but they are 
unequal in general. Though it does not seem easy to obtain the individual p„(w)'s for 
different Catalan words, it is, however possible to calculate their total contribution 
in a relatively simple way. We turn to this direction now. 
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Table 1. Pu{w) for Catalan words for 



Word 


Pu{w) 


dbdb 


1/2 


aabb 


1/3 


abba 


1/3 


aabbcc 


1/4 


abbcca 


1/4 


abbacc 


5/24 


aabccb 


5/24 


abccba 


5/24 



Theorem 3.1. Let be a triangular Wigner matrix with an input sequence satis 
fying Assumption C. Then almost surely 

lim — Tr 

n— s>oo Ti 



' ' n 



Hence the density of the LSD is given by where satisfies (1.5) and its 

support is contained in [—y/e, y/e\. 

Note that since the LSD is universal, the above limit moments are also the moments 
of the LSD under Assumption A. 

Remark 3.1. Interestingly, the above moment sequence is directly related to the Lam- 
bert W function, which is defined through the following functional equation 

W{z)exp{W{z)) = z. 

If X is real then the above equation has two possible branches and the one satisfying 
W[x) > —1 is called the principal branch and is denoted by Wq. The principal branch 
Wo is analytic at zero and one can use Lagrange's inversion theorem to see that, 

1 ^ k'^ 



1 + 



^ (k 

k=0 ^ 



-X 



xWo{l/x) 

for X lying outside (— e,e). For details we refer the reader to the exciting article on 
Lambert W function by Corless et al. |10j . 

Proof. To begin with, we need the following the simple lemma whose proof is omitted. 

Lemma 3.1. (i) Let w be a Catalan word of length 2k. Let n be a function vr : 
{0,1,. ^ {l,2,...,n} such thatw[i] = w[j] (7r(i - l),7r(i)) = (7r(j),7r(j - 1)). 
Then tt{0) = Ti{2k) and hence vr G 11* [w). 

(a) Let (f)(j) be as in Definition \2. 1 If w = W1W2 where Wi and W2 are Catalan 
words of length 2ki and 2k2 respectively. Then for every vertex i > 2ki, (f){i) > 2ki. 

Now let w he a Catalan word of length 2k. Let S be the set of generating vertices 



for w. It has already been shown that (see (2.8)) 

Pu{w) = ■■■ / I{v^(i^i) + v^(i) <l,ieS- {0})dvs. 



(3.1) 



[0,1] 



fc+i 



TRIANGULAR MATRICES 



13 



This integral can be evaluated as an iterated integral. It is clear that integrating out 
all the variables other that Vq leaves a polynomial in vq, say Q^- It follows that 

Pu{w) = / Qwivo)dvo (3.2) 







for some polynomial Qw{')- Let us illustrate with a few examples. Let w = aa. Then 

Pu{w) = / dvodvi= / {l-Vo)dvo 

Jvo+vi<l Jo 



and hence 

Now let w = abba. Then 



Hence 



J " J -^("^0 + "i^i < 1, "i^i + "^^2 < l)dvodvidv2 

[0,1]3 

/ / dv2dvidvo 
Jo Jo 

1 /•l-t'O /"I n _ y2\ 

{1 — vi)dvidvo = / ° (ifQ. 



1 - 



~'0 



2 

The following two lemmas collect the required properties of Q^(-). 

Lemma 3.2. (i) Let w = wiW2 be a Catalan word of length 2k where wi and W2 are 
both Catalan. Then Qw{x) = Qwi{x)Qw2{x) . 

(a) Let w = awia be a Catalan word with wi Catalan. Then 

pl—x 

Qw{x) = / Q^,{y)dy. 
Jo 

Proof, (i) Let wi and W2 be Catalan words of length 2ki and 2A;2 respectively. We 
divide the set of inequalities in the indicator function in (3.1) above into two classes: 

(1) + f < l,i E S — {0},i < 2ki; i.e., the inequalities corresponding to 
generating vertices in wi. 

(2) + f < l,i E S — {0},i > 2ki; i.e., the inequalities corresponding to 
generating vertices in W2. 

From Lemma [3.1| (i) it follows that (p{2ki) = and from Lemma 3.1 (ii) it follows 
that the inequalities in items (1) and (2) do not have any common variable except 
Vq. Hence in (3.1) integrating w.r.t the variables corresponding to the generating ver- 
tices in wi we get Qwi{vo) and integrating w.r.t. the variables corresponding to the 
generating vertices in W2, we get Qiu2(^o)- The result now follows from the definition 
of Qto- 



(ii) The proof of this is similar and we only give a sketch. Note that the variable vq 
does not occur anywhere in the equations corresponding to the generating vertices in 
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Wi. Integrating w.r.t. all these variables (other than t>i), we get Qwii^i). Hence the 
last step of evaluating the iterated integral is by ( |3.1 ) 



Pu{W) = 

and the result follows. 
Let 



vo+viKl Jo Uo 



dv, 



'0, 



Go{x) = 1 and G 



2n[X 



w Catalan, 
|ui|=?i 



X . 



□ 



(3.3) 



Lemma 3.3. With {G2n,'n > 0} as in (3.3), 

(i) for n > 1 we have 



G'2n(a;) = V'G'2(n-fe)(a;) / G2{k-i){y)dy. 

Jo 



k=l 



1-x 



(n) 



GoJx) 



[1 - x){n + 1 - x)""-^ 



for all n > 0. 

Proof, (i) For n = 1, the only Catalan word of length 2 is aa. Hence G2{x) = 
Qaa{x) = 1 — X and the result is true for n = 1. Let n > 2 and let G2n,k{x) be the 
sum of Qw{x) over all Catalan words w such that the first letter is repeated at the 
2/c-th place. Clearly, for such a. w, w = awiaw2 where wi is a Catalan word of length 
2{k — 1) and W2 is a Catalan word of length 2{n — k). Using the previous lemma, it 
follows that 



G2n.k(x) 



X 



\wi\ = {k-l),\w2\={n-k) 

^ ^ Q aw\a(jX:^Qw2 

\wi\ = (k-l),\w2\={n-k) 

\wi\ = (k-l) \w2\={n-k) 
pl—x 

/ QwAy)dy 



\wi\ = {k-l) 



\w2\ = {n~k) 



W2 



X 



QwAy))dy 



\w2\={n-k) 



ki|=(fe-l) 



G2(n~k){x) / G2(k-i){y)dy. 
Jo 



As G2n{,x) = Ylk=iG2n,k{x), part (i) follows. 
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(ii) We prove this by induction. The cases n = 0, 1 are clear. Now suppose that the 
result is true for all j < n. Then by the last lemma and the induction hypothesis we 
have 

-l-x 



pi—x 

G2n{x) = y'G'2(n-fc)(a;) / ^2(^-1) 

^0 



k=l 

n 



^ (1 - x){n -k + 1- a;)"-^'-i f^''' (1 - y){k - yf'"^ 



k=l 



-^{l-x){n-k + l- x)"-''-^ (1 - a;) (A; - 1 + xf'^ 



k=l 



{n-k)\ k\ 

(^Pk{z)Pn-k{-z), 



k=l 

where z = 1 — x and 

Pn{x) = x{x — n)^~^ 

is the Abel Polynomial of degree n. It is a well known fact (see Riordan [16]) that 
Abel polynomials satisfy the following combinatorial identity: 



Pn{x + y) = Y^ r^jPk{x)Pn-k{y)- 

k=0 ^ ^ 



It follows that 



G2n{x) = j 22 [ )^]Pkiz)Pn-k{-z) + ^Po{z)Pn{-z) 



k=0 

(-1)" 

nl nl 
(-l)"(-2)(-2 - 



z{z + n) 



n-l 



{I - x){n + I - xy-^ 
n\ 

and the proof of Lemma 3.3| is complete. □ 
Now we complete the proof of Theorem |3.1[ Using the previous lemmas we have 

lim -Tr(— ^) = > Pu{w) 

w Catalan, 
|ui|=fc 

1 



G2k{x)dx (see equation(3.3)) 
^ {1 - x){k + 1 - x)''-^ 



k\ 



dx (Lemma 3.3K ii)) 
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□ 



4. Some comments on other variants and joint convergence 

(1) It can be shown that the hmit which appears in the triangular Wigner matrix 
can appear in a variety of other matrices which satisfy the so called Property 
P introduced in [2] . The Property P states that 

M* = sup sup |{1 < A; < n : + i < n + 1, A; + j < n + 1, L{k,i) = L{kJ)}\ < oo. 



The proof of this fact is similar to the proofs in in [2] and hence we skip it. 

(2) The joint convergence of patterned random matrices was initiated in [Hj and 
was further studied by [5] for multiple copies of different independent pat- 
terned random matrices. We can so study the joint convergence of indepen- 
dent patterned triangular matrices along the same lines. It can be shown that 
if the entries have uniformly bounded moments and if Puiw) exists then in- 
dependent copies of patterned random matrices jointly converge with respect 
to - Tr for any of the matrices considered in this article. The behavior in the 
limit though, is a highly non-trivial problem. 

(3) Let be two triangular Wigner matrices having uniformly bounded 
moments. Let ({ai,a2},0) denote the joint limit of ({^7=^; }; 0n) where 
<^n = ^E[Tr(-)]. Then ai and 02 are not free. In fact, 



Urn -E[TT{W^jV^y] =pu,w{aa) = 1/2. 

1 — Vrvn T) ' 



Also, from Table 1, 




(4) Another variant of the triangular Wigner matrix is 






















X2{n-l) 



w 



n 








•^(n— l)(n— 1) -^(n— l)n 
•^{n—l)n •^nn 



n 



It should be noted that we can always delete the diagonal from the above 
matrix without changing the limiting spectral distribution. Let P be the 
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following orthogonal matrix 



Then we have 



pw::p' 














1 








1 









X2(n-1) 





1 








Xin 
Xi2 Xii 



(4.1) 



(4.2) 



which can be seen to be equal to upon renaming of the variables. Since 
conjugation by orthogonal matrices preserves eigenvalues, we can conclude 
that the limiting spectral distribution of W}^/ y/n is same as that of W^. 

(5) Interestingly, W^/y/n and VF^/i/n jointly converge. If we call the hmiting 
variables as 01,02, then again it can be shown that they are not free. But it 
easily follows that Oi + 02 obeys the the semicircular law. 



(6) As is well known, for the full version of the matrices, the contribution of every 
Catalan word to the limiting moment equals one for all the common patterned 
matrices. Further, for the Symmetric Circulant every word contributes one. 
For the Reverse Circulant, only the so called symmetric words contribute, each 
contributing one. For the Hankel matrices also, only the symmetric words con- 
tribute but all of them do not contribute one or even contribute equally. See 
for example |8]. It is an interesting question to study the individual contri- 
bution of each word to the limiting moment for the triangular matrices. This 
does not appear to be easy, even for the Wigner case. 

(7) The following facts on word contribution can be proved for the triangular 
Wigner matrices. The proof is given in Section [5j 

(i) For any word w of length 2k, Pu{w) < l/{k + 1). 

(ii) If w is of the form w = aabbcc. . . then Pu{w) = l/{k + 1). 

(iii) If w and w' are Catalan words of the form w = awiaw2 and w' = aawiW2 
then Pu{w) < Pu{w'). 

(iv) If w and w' are Catalan words of the form w = ahhawiW2 and w' = 
abwibaw2 are Catalan words then Pu{w) > Pu{w'). 

(8) We wish to point out some histogram plots for the patterned triangular ma- 
trices. In Figures |2] and [3] we have plotted the histograms of the empirical 
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spectral distribution of the symmetric triangular Wigner, Toeplitz, Hankel 
and Symmetric Circulant matrices. 



It can be easily checked from the density equation given earlier in |1.5| that 
the density is unbounded at zero for the triangular Wigner. Further, from the 
moment formula, one can easily check that the support of the LSD is contained 
in [— i/e, \fe\. We can see evidence of both these facts in the histogram. It is 
known that the maximimum eigenvalue of the full Wigner matrix converges 
almost surely to 2. We believe that likewise, for the triangular Wigner, the 
maximum eigenvalue converges to -y/e almost surely. 

As pointed out earlier, obtaining moment properties of the LSD for the 
other three matrices does not seem to be easy. However, it appears from the 
simulations that all the other three LSDs also have densities and these are 
also unbounded at zero. Moreover all these LSDs have unbounded support, 
just like for their full matrix versions. 

Interestingly, it is known that the LSD of the full Hankel is not unimodal 
(simulations results show that it is bimodal but there is no formal proof). 
However, the simulation evidence implies that the triangular Hankel LSD is 
unimodal. Recall that for the full Hankel matrix, the entries in the upper 
triangle and the entries in the lower triangle (ignoring the main antidiagonal 
part) are independent of each other. One wonders whether this particular 
nature has anything to do with the bimodality of the LSD of the full Hankel. 



(9) Now suppose we consider triangular patterned matrices but we drop the as- 
sumption of symmetry. Let be any such matrix. It can be shown by the 
moment method that the LSD of the symmetric matrix A^A'^/n exists. A 
natural (and hard) question is what happens to the LSD for the asymmetric 
matrix Anj \/n. The moment method cannot be directly applied and one has 
to resort to the Stieltjes transform. 

(10) Let X^^^ be a symmetric matrix with link function L and let X be it asym- 
metric version as defined in [6]. A general question that was raised there was 
the following. Suppose X^^^ / ^Jn has an LSD identified by a random variable 
say X. Suppose XX' jn has an LSD denoted by Y say. When do and Y 
have the same distribution? 

The answer is affirmative when X^^^ is a symmetric Wigner matrix and 
X is its asymmetric version (the matrix with i.i.d. entries) from the relation 
between Marcenko-Pastur law and the Wigner law. In [Sj it was shown that 
this is true for the Toeplitz matrix but not for Hankel or Reverse Circulant 
matrices. 

The result in this article gives another important example in the form of 
the triangular Wigner matrix. 
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Figure 2. Histo gram plots of empirical distribution of triangular matrices (n = 
1000) with entries N{0, 1) of (i) Wigner (left) (ii) Toeplitz (right) 




Figure 3. Histo gram plots of empirical distribution of triangular matrices (n = 
1000) with entries N{0, 1) of (i) Symmetric Circulant (left) (ii) Hankel (right) 

5. Appendix 

5.1. Contribution of different Catalan Words for triangular Wigner Matrix. 

It appears difficult, even for triangular Wigner matrix, to determine Pu{w) for each 
Catalan word w. However, Pu{w) can be determined for some special classes of words. 
As we have already seen, unlike what happens for a full Wigner matrix, Pu{w) is not 
same for all Catalan words. Here we record some observations about contributions 
of different Catalan words for triangular Wigner matrices. This might be useful in 
obtaining more information about the individual j9„(w)'s. We start with a simple 
lemma. 
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Lemma 5.1. Let w be the Catalan word of length 2k, w — aiaia2a2...akak. Then 

Proof. Clearly, tt G 11* (w) implies 

(min(7r(0), 7r(l)), max(7r(0), 7r(l))) = (min(7r(l), 7r(2)), max(7r(l), 7r(2))) 
and hence 7r(0) = 7r(2). Arguing similarly, it follows that. 



n*{w) = {n : 7r(0) = n{2) = 7r(4) 



7l{2k)}. 



Hence, 



|n*(w) I = #{7r : 7r(0) + 7r(l) < n, n{0) + tt{2) < n, . . . , , n{0) + n{k) < n}. 



Now 



lim 



\ni{w) 



n— >oo 77, 



1+fe 



J "j + ^"1 < 1, ^^0 + ^"2 < 2, Vo + Vs< 1, Vo + Vk< l)dvodvi...dvk 

nl r rl-vo 

Jo Uo 



dvi ■ ■ ■ dvk 



VQ^dvo 



1 



k + l 

which completes the proof. □ 

Now we show that words in some classes contribute more than words in some other 
classes. Before that wc need the following two lemmas. 

Lemma 5.2. Let w he a Catalan word of length 2k. Let tt be a function tt : 
{0, l,...,/c} ^ {l,2,...,n} such that w\i] = w\j] =^ (7r(i - l),7r(i)) = (7r(j),7r(j - 1)). 
Then 7r(0) = 7i{2k) and hence n G Il*{w). 



Proof. The proof is easy by induction and we omit the details. 



□ 



Lemma 5.3. Let w be a Catalan word of length 2k such that w = awiaw2, where Wi 
and W2 are Catalan words. Let w' denote the Catalan word aawiW2. Then Pu{w) < 
Pu{w'). 

Proof. Let us denote, for any set U G M'^'^"'^, by Vol{U), the Lebesgue measure of the 
set [^nP) 1]^^^- From what we have already proved, it follows that for every Catalan 
word wq of length 2k, Pu{wo) = VoI{Uq) where Uq is the region in R*^"*"^ given by the 
set of all v's determined by the set of inequalities 

N = K(i-i) + v^^i) <l,ieS- {0}}. (5.1) 

Now let w — awiaw2 and let wi and W2 be Catalan words of lengths 2ki and 2k2. 
Now we partition the set N into three classes. 
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(1) Ni = {vq + vi < 1}: This is the first inequahty since 0(0) = 0,0(1) = 1 and 
1 e S. 

(2) Aw-j^{vi): this is the set of all inequalities + f<^(j) < I ■ 2 < i < 
2ki + l,i G S — {0}}; i.e. the inequalities corresponding to the generating 
vertices in wi. By Lemma [3. 1| it follows that these inequalities do not involve 
the variable vq. Also, the inequalities only depend on wi and vi. 

(3) Aw^{vq): this is the set of all the inequalities {f<)<,(j_i) +f</,(i) < 1 : 2ki + 3 < i < 
2k, i G S* — {0}}; i.e. the class of inequalities corresponding to the generating 



vertices in W2- From Lemma 5.2, it follows that 7r{2k + 2) = 7r(0) for every 



TT G Il*{w) and hence (f){2k + 2) =0. Now Lemma pA implies that all the 
variables occurring in these inequalities are independent of the variables that 
occurred previously, apart from t>o. Also, the inequalities only depend on wi 
and Vq. 

It follows that Pu{w) = Vol{Bw) where 

= NiUA^^{vi)UA^^{vo). 
Similarly, for the word w' = aawiW2 the three classes of inequalities are: 

(1) N,: 

(2) Aw-^{vq): in this case 0(2) = 0. By renaming the variables other than vq, if 
necessary, the class of inequalities corresponding to the generating vertices in 
Wi in this case is same as that in the previous case with vq replaced by vi. By 
Lemma |3.1 it follows that these inequalities do not involve the variable Vi . 



(3) Ay^^{vo): the class of inequalities corresponding to the generating vertices in W2 
in this case is same as that in the previous case after renaming the variables 



other than vq- Lemma 3A_ implies that all the variables occurring in these 
inequalities are independent of the variables occurred previously, apart from 

Vq. 

It follows that Pu{w) = Vol{Byj') where 

B^> = NiUA^^{vo)UA^^{vo). 



Now 



where 



Vol{B^) = Vol{Bl) + Vol{Bl) and 
Vol{B^.) = Vol{Bl,) + Vol{Bl,) 



Bl, = NiUA^,{vi)UA^,{vo)U{vo<vi}, 
Bl = NiU A^,{vi) U A^,{vo) U {vi < vo} 



and 



Bl, = NiUA^,{vo)UA^,{vo)U{vo<vi}, 
Bl, = N^U A^^ (vo) U A^^ (vo) U {v^ < vq}. 

It is easy to see that Bl, D Bl and Bl D Bl,. 
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Now, 

Bl, -Bl = iVi U A^^{vo) U A^^{vo) U{vo<v,}U (A^.K))'; 

Bl-Bl, = niUA^,{vi)UA^,{vo)U{vi<vo}UiA,,„{voy. 

Now by interchange of variables Vq and Vi, it follows that Vol{B^, — B^) = Vol{Cu,,w') 
where 

C^,«;' = A^i U A^^{vi) U ^^^(^^i) U < vo} U (^^,(t;o))'. 
Again, vi < vo and Atu2(i'o) together imply Auj^ivi) < 1 and hence Cu,,w' ^ B^j—B^,. 
It follows that, 

l^oi(C^,«,') > Vol{Bl-Bl,) 
^Vol{Bl,-Bi) > Vol(Bl-Bl,) 
^Vol{Bi,) + Vol{Bl) > Vol{Bl) + VolBl) 

and this completes the proof. 

□ 

The next Lemma is similar to the previous one. 

Lemma 5.4. Let w = abbawiW2 and w' = abwibaw2 be two Catalan words where wi 
and W2 are Catalan words. Then pu{w) >Pu{w'). 

Proof. We use the same notations. As in the previous Lemma, now Pu{w) = V^o/(i?^) 
and Puiw') = Vol{Byj') where 

By, = A^i U + < 1} U A^^ (vq) U A^^ {vo) and 
By,, = Ni\j{vi + V2}UAy,^{v2)UAy,^{vo). 

Now 

VoliBy,) = VoliBl) + Vol{Bl) and 
Vol{By,,) = Vol{Bl,) + Vol{Bl,) 

where 

Bl = NiU {vi + V2<1}U A^,, (i'o) U A,,, (t-o) U {^;o < V2}, 
Bl = NiU {vi + V2<1}U Ay,^ {vo) U Ay,^ {vo) U {v2 < Vo}, 

and 

Bl, ^ NiU {vi + V2<1}U Ay,, {V2) U Ay,, (vo) U {vo < V2}, 
Bl, ^ NiU {vi + V2<1}U Ay,, {V2) U ^1^2(^0) U {V2 < Vo}. 

As Vq < V2 and Ay„{v2) together imply Ay,,{vo), Bl D Bl,. Similarly it follows 
that Bl, D Bl. 
Now, 

Bl - Bl, ^ NiU {vi + < 1} U Ay„{vo) U Ay„{vo){vo < V2} U {Ay,^)'; 
Bl, -Bl = NiU {vi + < 1} U Ay„{v2) U Ay„{vo){v2 < Vo} U {Ay„{vo)y. 
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Once again by interchange of two variables Vq and V2 we see that Vol{Bl^, — 5^) = 
Vol{Cii,,w') where 

C^X = A^i U + ?;2 < 1} U A^,{vo) U A^,{v2){vo < V2} U {A^,{v2y . 
Again, vq < V2 and ^4^2(^2) together imply 74^^,2(^0) and hence Cw,w' C 5^ — B^,. 
It follows that, 

Vol{Cw,w') < Vol{Bl-Bl,) 
^Vol{Bl,-Bl) > Vol{Bl-Bl,) 

and the proof is completed as in the previous Lemma. □ 

The next proposition gives an upper bound on pu^wYs. 
Proposition 5.1. Let w be a Catalan word of length 2k. Then Pu{w) < 



Proof, li w = 01010202. ..afcafc, then the result follows from Lemma 5.1 If not, by left 
rotation, we can obtain, from w, a word w' such that w' does not start with a double 
letter and w' is not of the form owio either. Since Pu{w) is invariant under rotation, 
Pu{w) = Puiw'). By hypothesis, w' must be of the form w' = awiaw2 where Wi and 



W2 are Catalan words. By Lemma 5.3, Pu{w) = Pu{w') < pu{aawiW2) = Pu{wiW2aa) 
and the number of consecutive double letters at the end of ^1^200 is strictly greater 
than that of w. The proof can now be completed by induction on the number of 
consecutive double letters at the end of w. □ 

We can see from the above results that Puiw) depends upon the Noncrossing struc- 
ture of the word w. It is an interesting combinatorial problem to obtain a formula 
for Pu{w) as a function of the Catalan structure of w or to investigate what Catalan 
structures gives rise to the same pu(w)'s. 
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